Enhanced radiative heat transfer between nanostructured gold plates 
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We compute the radiative heat transfer between nanostructured gold plates in the framework of 
the scattering theory. We predict an enhancement of the heat transfer as we increase the depth of 
the corrugations while keeping the distance of closest approach fixed. We interpret this effect in 
terms of the evolution of plasmonic and guided modes as a function of the grating's geometry. 



PACS numbers: 64.70.Nd, 44.40. +a, 44.05.+e 

The far-field radiative heat transfer between good con- 
ductive metals is very low at room temperature, since 
they are very good reflectors at the infrared frequencies 
of blackbody radiation. The radiative heat transfer is 
enhanced in the near field, due to the contribution of 
evanescent surface modes |IH3] ■ Polar materials like Si02 
or SiC are in addition favored by the contribution of sur- 
face phonon polaritons whose resonance frequencies lie in 
the infrared [1]. There is an analogous effect for metals 
arising from the surface plasmons resonances but those 
lie in the ultraviolet and do not contribute significantly 
to the heat transfer [S]. 

It has been shown recently that the radiative heat 
transfer can be controlled by nanostructuring the inter- 
faces periodically. When the period d is much smaller 
than the wavelength A and the separation distance L, 
the system can be treated using an effective refractive in- 
dex for the equivalent homogeneous medium. It has been 
shown that the induced anisotropy introduces additional 
modes [B] and also allows modulating the fiux [7]. For 
periods on the order of the wavelength, a full solution of 
Maxwell equations is needed. The heat transfer between 
two periodic slabs has been studied within a two dimen- 
sional approximation and for p-polarization using a finite 
difference time domain (FDTD) technique [5]. A fiux en- 
hancement attributed to the excitation of the structure's 
modes was found. While FDTD allows modeling complex 
shapes easily, dealing with bulk 3D media and accounting 
for polarization effects has not been achieved so far. 

In this letter, we compute the radiative heat trans- 
fer between ID gold lamellar gratings in the framework 
of the scattering theory. We do include all propagation 
directions (the so-called conical diffraction) and all po- 
larization states, which is of critical importance in order 
to deal quantitatively with cross-polarization effects [3]. 
The scattering theory is the most successful technique for 
treating the Casimir effect between bodies at thermody- 
namic equilibrium |101 lllj . The method determines the 
electromagnetic field in the space between the two bod- 
ies in interaction in order to compute the Maxwell stress 



tensor in terms of the reflection amplitudes on the two 
bodies. When the two bodies are not at the same temper- 
ature, there is a net fiux of energy transferred from the 
warm body to the cold one. Recently, this heat transfer 
problem between two bodies kept at different tempera- 
tures has also been formulated in terms of the scattering 
properties of the bodies P^HTS] . 

In the following, we use the scattering amplitudes 
which have already been calculated for studying the 
Casimir interaction between ID lamellar gratings [TS] 
and deduce the heat fiux when the two bodies are at 
different temperatures. We show that the heat flux is 
largely enhanced when the corrugation depth is increased 
while keeping the distance of closest approach fixed. We 
attribute the heat fiux increase to the excitation of guided 
modes and surface plasmons whose frequencies change 
with the corrugation depth. 
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FIG. 1. The conventions used in the present paper. The 
grating period is d, the corrugation depth is a, the distance 
of closest approach of the two gratings is L. The lines of the 
grating are along the y direction, while the Fabry-Perot cavity 
between the two gratings is along the z direction. 

We consider the cavity formed by two gratings sepa- 
rated by a distance of closest approach L measured so as 
to vanish at contact (Fig. [I]) . The gratings are aligned 



and not displaced laterally. We model the gold permit- 

2 

tivity with a Drude model e(w) = 1 7-^ — r with wp = 9 

eV and 7 = 35 meV. We write the heat flux q between 
two bodies at temperatures Ti and T2 as [3 



q 



(eTiH-eT2(w))rL(k,a;)— -p-, (1) 



where eT{<-ij) — fujj [gf^l^sT _ ^j jg ^j^g niean energy 
per mode of frequency w at temperature T while TKk, lS) 
is the sum (trace) of the transmission factors for all the 
modes of frequency w and lateral wavevector k between 
the two gratings separated by a distance L [T71 [THj . The 
expression of this transmission factor is given by scatter- 
ing amplitudes 

Tl (k, u^) = tr (DWiDt W2) , (2a) 

D = (l-SiS2)-\ (2b) 

Wi = S^"; - SiS^^Si^ + SiSl'"! - Sl'^iSi^ (2c) 

W2 = Sr - S2^SrS2 + S2^Sr - S^Sa, (2d) 
Si=Ri(k,w), (2e) 
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Mode counting is defined over frequency uj and lateral 
wavevector k restricted to the first Brillouin zone, due to 
the Bloch theorem, k^ = ^llP- jc? — k^ is the longitudi- 
nal wavevector for the Fabry-Perot cavity, with the prin- 
cipal square root used in its definition — f < arg A:^ < ^ . 

The operators H^''^^ — J^^YIP'^/^"' involve the projec- 
tors nP™/^™ on the propagative or the evanescent sector, 
respectively. Si and S2 are scattering operators defined 
from the reflection operators Ri(k,a;) and R2(k, w). Si 
are represented in the basis of the wavevectors {k^")} 
coupled by the grating. We define k*^"^ = k -I- «.^e^ 
where d is the grating period, e^, the direction perpen- 
dicular to the lines of the grating (see figurefl]) and n runs 
from —N to -l-A^, where N is the highest diffraction order 
retained. The operators Si are square matrices of dimen- 
sion 2 (2 A^ -1-1) [12] as well as all bold operators appearing 
in eqs.[2] All scattering operators appearing in eqs.[2]are 
represented in the (s/p) (also denoted TE/TM) polar- 
ization basis, well adapted to propagative fields. The 
reflection operators are calculated following the Rigor- 
ous Coupled- Wave Analysis (RCWA) method described 
in [12]: the fields are expressed in terms of a Rayleigh 
expansion in both homogeneous regions z < and z > a. 
In the corrugated region > z > a, the fields are de- 
veloped in Fourier components. The Maxwell equations 
are solved in each regions and writing the continuity of 
each Rayleigh and Fourier components at the boundaries 
z = and z — a leads to the reflection and transmission 
coefficients for the grating. In the limit of an infinite 
number of Fourier harmonics, this method solves exactly 
the diffraction of the fields by the grating. Metallic grat- 
ings are known to be difficult to account for using the 



RCWA method. We incorporate in the RCWA formal- 
ism the modifications presented in |20j which greatly im- 
prove the convergence rate for the reflection coefficients 
of a p-polarized light impinging on a metallic grating and 
our calculation are performed with A'^ = 51 which shows 
converged results. 

In the following, we apply formula (fTj) to compute the 
heat transfer coefficient h defined as ft. = rp ^^ for two 
temperatures Ti and T2 close enough to each other, say 
for example Ti = 310 K and T2 = 290 K. We note that 
gti ~ eT2 9'Cts as a cutoff function for frequencies greater 
than the thermal frequency lot — x^ ~ 2.5 x 10^'' rad s~^ 
(At ~ 7.6 /im). The transmission factor 7L(k,a;) thus 
exhibits the mode structure for the problem under study 
(Fig. [1]) while ([l]) integrates the contributions of all these 
modes to the heat transfer, taking into account the values 
of their frequencies with respect to lot (more discussions 
below). 

For a depth of the corrugation a = 0, we recover the 
heat transfer coefhcient h^^L) — 0.16 Wm~^K~^ be- 
tween two gold plates separated by a distance L = 1 
/im. For a non null depth a, we introduce the factor of 
enhancement of heat transfer with respect to non corru- 
gated plates 
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h{L) 
ho{L)- 



(3) 



We present in Fig. [2] the enhancement factor fi as a func- 
tion of the corrugation depth a, with the distance of clos- 
est approach L = 1 /xm and the filling factor p = 0.5 kept 
fixed. The blue solid curve corresponds to a period d = 1 
/im for the gratings while the red solid curve corresponds 
to a period d — 2.5 /im. The dashed curve corresponds 
to a period d = 10 /im. As the corrugations become 
deeper, we see a striking increase in the heat transfer co- 
efficient. We note that the enhancement factor is largely 
independent of the grating period up to a corrugation 
depth awl /xm. For a period d = 1 /xm for which the 
effect is more important, we get an enhancement up to 
a factor 10 for a ~ 6 /im. For a period d = 2.5 /im, 
the enhancement reaches nearly a factor 4 for a = 6 /im. 
For the largest period d = 10 /im, the enhancement still 
reaches nearly a factor 2 at a = 6 /im. 

For comparison, we have shown as the dotted line 
in Fig. [2] the prediction of the proximity approximation 
(PA) which amounts to adding plane-plane heat transfer 
contributions, as if they were independent. 
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PA 



P+{l-p) 



ho{L + 2a) 
ho{L) 



(4) 



As expected, the PA predicts a decrease of H. when a 
is increased, in complete contradiction with the exact 
results shown by the solid and dashed curves. 

In the remainder of this letter, we analyze the electro- 
magnetic mode structure in order to explain the increase 
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FIG. 2. The enhancement factor Q, between two gold gratings 
as a function of the depth a of the corrugations, with the 
distance of closest approach kept fixed L = \ fim. Red solid 
curve: period d — 2.5 fim. Green dashed curve: period d = 10 
fj,m. Black dotted curve: proximity approximation, (colors 
online) 




of the heat transfer [T71 HH]. To this aim, we use the 
scattering formula (IT]) and show that, as we increase the 
corrugation depth, some modes of the system are indeed 
brought to the infrared frequencies and thus are able to 
contribute to the heat transfer. The mode structure is 
described by the transmission factor 71 (k, w) which can 
reach its maximum value 1 at the resonances of the cor- 
rugated cavity. Our system is periodic so that the mode 
structure, distributed over the whole range of wavevec- 
tors in the absence of corrugations, now shows many 
branches folded in the first Brillouin zone. More pre- 
cisely, there are 2(2N + 1) branches where the factor 2 is 
due to the two polarizations and the factor 2N -I- 1 is the 
number of orders (or branches) used when taking into 
account mode coupling by diffraction on the gratings. 

We represent in Fig. [3] the sum of transmission factors 
TiCk, ijj) over all polarizations and all branches. It is 
shown as a function of the frequency uj and the depth 
of the corrugations a for a fixed value of the transverse 
wavevector k = (53, 0), here chosen to be in the middle of 
the positive-fca; first Brillouin zone. The plot corresponds 
to the period d = 2.5 /xm, which was shown as the solid 
curve in Fig. [21 The vertical red line represents the light 
line oj = ckx « 1.88 x 10^** rads^^. 

It clearly appears in Fig. [3] that the transmission factor 
takes significant values only on resonances which corre- 
spond to the mode structure of the corrugated cavity. 
The transmission factor 7l (k, w) goes to a maximum 
value of 1 for each non degenerate mode (k, oj); it can 
be 2 if two modes cross and we see one of these occur- 
rences in the figure. The general trend is clear on the 
diagram: as the depth a of the corrugations is increased, 
new modes appear, with frequencies decreasing as a in- 
creases. When these modes enter into the thermal win- 
dow UJ < ujt they contribute more and more to the heat 



FIG. 3. The transmission factor for two gold gratings as a 
function of the frequency uu and the corrugations depth a. 
The lower curve is for plane-plane a — while the upper one 
is for a corrugations depth a = 3 /im. The vertical red line is 
the light line. The horizontal arrow at a = 1.5 fim shows a 
cut of this plot represented on Fig. B] (colors online) 



transfer. This explains the enhancement of the heat flux, 
due to the presence of additional modes in the thermal 
window for a deeply corrugated structure. 

We now examine in more detail the nature of the 
modes. While varying the corrugation depth a from 
to 3 /im we can follow the evolution of each mode. Note 
that, for ky = 0, the polarizations a = s and a = p 
are not mixed (however, the computation of h takes into 
account all modes for which polarization mixing is im- 
portant). 

We show in Fig. |4] the modes calculated for a partic- 
ular corrugation depth a = 1.5 /im indicated by the red 
horizontal line on Fig. [3] The position of the peaks have 
been confirmed through a direct mode calculation |21j 
of the eigenfrequencies of the structure modes obtained 
for p (black arrows) and s (red arrows) polarizations. In 
addition to the excellent agreement between the peaks of 
the transmission factor and the directly calculated modes 
(arrows on Fig. Ill), direct mode calculations show the 
fields and, therefore, allow us identifying the first few 
modes. For the second p polarization and the first s po- 
larization modes appearing at w « 2.4 x lO^'* rads~^ and 
Lu K, 6.5 X 10^^ rads~^ in particular, the frequencies are 
largely independent upon the value of fcx, which is usu- 
ally the signature of guided modes. By looking at the 
fields corresponding to those two modes, we indeed con- 
firmed that the electric field is to some extent confined 
in the waveguides formed by the corrugations. 
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FIG. 4. The transmission factor for two gold gratings with 
corrugation depth a — 1.5 fim as a function of frequency 
CO. The arrows indicate the position of the modes in a direct 
mode calculation (red for s polarization and black for p polar- 
ization). The dashed curve is the function ^ Yt -t ) ■ (colors 
online) 
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FIG. 5. Same as fig. [4] for the first two modes which are in 
the thermal window. Each peak can be fitted by a Lorentzian 
of resonance frequency wq and half width at half maximum 
F. The dashed curve is the function -. — tts — ^. 



considered here of lossy materials. 

This discussion allows one predicting the effect of a 
change of the dissipation parameter 7. As this parameter 
is the only one to determine the widths of the peaks in 
the transmission factor TlO^,oj), one deduces that these 
widths vary hnearly with 7. As a direct consequence of 
(IT]) and as long as the modes remain sharp and isolated, 
it follows that the heat fluxes vary in proportion of 7, so 
that the enhancement factor fJ, deflned in (Is]) and drawn 
on Fig. [2] is independent of the dissipation parameter 7. 

We have theoretically demonstrated the enhancement 
of the heat transfer between two nanocorrugated gold 
plates in comparison with flat plates with the same dis- 
tance of closest approach. This enhancement is due to the 
presence of additional modes in the thermal frequency 
window contributing to the heat transfer. We have 
described all the relevant information about the mode 
structure in terms of the transmission factor 7i(k,w) 
which appears in the scattering formula for the heat flux. 
We have discussed the enhancement of the heat transfer 
in a regime where the three characteristic lengths of the 
problem (the distance L between the gratings, the period 
d of the gratings and the height a of the corrugations) 
are of the same order. We stress that neither the prox- 
imity nor the effective medium approximations can work 
in this regime. We have in fact shown that the proximity 
approximation predicts a decrease of the heat transfer, 
in complete contradiction with the striking enhancement 
of the heat flux observed in the exact results. 
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It is also worth discussing the shape of the resonance 
curve drawn by the variation of the transmission factor in 
the vicinity of a mode. In fig. [5j we focus on the modes 
which lie inside the thermal window. In the case con- 
sidered here of sharp, isolated modes, the resonance of 
the transmission factor shows a Lorentzian profile. We 
have checked that the two parameters of this profile are 
identified respectively to the real and imaginary parts of 
the complex frequency, with mode calculation of the dis- 
sipative structure defined with complex frequencies and 
real wavevectors [3]. This proves that the variation of 
the transmission factor contains all the relevant informa- 
tion about the mode structure. Not only the frequencies 
but also their finite lifetime are well described in the case 
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